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Abstract 

The thermodynamic nature of two-dimensional vortex matter is studied through a du- 
ality analysis of the XY model over the square lattice with low uniform frustration. A 
phase-coherent vortex lattice state is found at low temperature if rigid translations are pro- 
hibited. It shows a non-zero phase rigidity that is degraded exclusively by the creation of 
dislocation pairs. The unbinding of such pairs causes the vortex lattice to simultaneously 
lose phase coherence and to melt at a continuous (Kosterlitz-Thouless) phase transition. 
General phase auto-correlation functions are also computed, and these results are used 
to argue for the existence of a continuous melting transition of vortex matter in layered 
superconductors. 

PACS Indices: 74.60.Ec, 61.72.Bd, 74.60.Ge, 61.72.Lk 
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The nature of the mixed phase in layered superconductors remains badly understood 
theoretically. 1 For example, the rich phase diagram shown by high-temperature supercon- 
ductors in external magnetic field continues to yield surprises. 2 These systems are layered 
and extremely type-II. The minimum theoretical description of a layered type-II super- 
conductor in external magnetic field is one that neglects both magnetic screening and 
Josephson coupling, which is equivalent to a stack of isolated two-dimensional (2D) XY 
models with uniform frustration. 3-5 A 2D vortex lattice state is believed to exist at low 
temperature in such case. 1 Indeed, the elastic medium description predicts that the 2D 
vortex lattice melts continuously into an intermediate hexatic phase that retains the six- 
fold rotational symmetry shown by the solid state. 6-8 On the contrary, calculations of the 
long-range phase correlations in 2D and three-dimensional (3D) vortex lattices that also 
employ the elastic medium approximation find no evidence for macroscopic phase coher- 
ence due to the appearance of an infrared divergence. 9 Controversy therefore surrounds 
the theoretical existence of a vortex lattice state in pure type-II superconductors. 10 ' 11 

In this Letter, we show that a phase- coherent vortex-lattice state is in fact possible at 
low temperature in two dimensions if rigid translations of the vortex lattice are prohibited 
by surface barriers. A direct calculation of the phase correlations in the 2D XY model 
with low uniform frustration demonstrates that the infrared divergence alluded to above 9 
can be removed in a natural and unambiguous way through an asymptotic regularization 
scheme. This result is confirmed independently by a direct calculation of the phase rigid- 
ity. The Villain approximation is employed throughout. 12-14 We also demonstrate that 
phase coherence is lost at the continuous melting transition as a result of the unbinding 
of dislocation pairs and in a manner that is parallel to the loss of the shear rigidity. 6 No- 
tably, the renormalization group associated with the phase rigidity in the vicinity of the 
2D melting transition is identical to that of the shear modulus. Last, we show how the 
difference between the intermediate hexatic phase and the vortex liquid phase is reflected 
by a subtle difference in the phase factor of the phase auto-correlation functions. 

The layered XY model with low uniform frustration provides a qualitatively cor- 
rect theoretical description of vortex matter deep inside of the mixed phase in extremely 
type-II layered superconductors, where magnetic screening can be neglected. 5 In the ab- 
sence of Josephson coupling, its thermodynamics is determined by the energy functional 
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E ( X y = —JJ2rv cos[0(r + jia) — 4>(r) — A^r)] for the superfluid kinetic energy of an iso- 
lated layer in terms of the superconducting phase </>(r), with local phase rigidity J and 
with vector potential A = (0,2/r/x/a). Here / is the density of vortices and a denotes 
the square lattice constant. Monte Carlo simulations of the above uniformly frustrated 
XY model indicate that melting of the low-temperature 2D vortex lattice occurs at a 
temperature /csT m = J/20 for low density of vortices, 3 /. Monte Carlo simulations of 
the corresponding non-neutral lattice Coulomb gas obtain a quantitatively similar phase 
diagram. 4 This indicates that the Villain approximation for the uniformly frustrated XY 
model in two dimensions, 12-14 which reduces to such a Coulomb gas ensemble, is valid at 
temperatures even as high or higher than the vortex-lattice melting transition. 

The Villain approximation itself for the XY model over the square lattice amounts 
to the expression (exp[z J2^p(r)(j)(f)]) = Z v \p]/Z v [0] for any generalized phase auto- 
correlation function set by an integer source field, p(r), in terms of the quotient of the 
corresponding partition functions 

{n(r)} V r r / 

Above, n M (r) is an integer field that ranges over the links, (r, /u), of the square lattice, 
while (3 = J/ksT. The regime of validity of the Villain approximation for the XY model 
is generally at low temperatures. 14 The conservation equation expressed by the 5- function 
factors above can be treated exactly. After a series of manipulations, one is led to a 
factorized expression for general phase auto-correlation functions in terms of independent 
spin- wave (sw) and vortex (vx) contributions: 13 ' 14 

( e *5>*) = C^W^^*') e- , /p X<fi '. (2) 

The spin- wave factor above is equal to an exponential C sw [p] = exp[— ?? sw 2) p(l)27rG( 2 )(l, 2)p(2)]| 
of the Greens function = —V -2 over the square lattice, and it decays algebraicly with 
a correlation exponent ?? sw = (27T/3)" 1 . Note that 13 ' 14 27rG( 2 )(l, 2) = \n(r' /r 12 ), with 
r' ~ a. The vortex contribution above represents an average over (integer) vorticity Q(R) 
on the dual lattice of points R. These vortices form a Coulomb gas ensemble with a uni- 
form background at a density /. The ensemble is weighted by the Gibbs distribution set 
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by the Coulomb energy 



£ vx = 27rJ^[Q(a)-/] 27vG( 2 \a,b) [Q(b) - /]. (3) 

(a,6) 

The flux field 

$ P (R) = -2nJ2 G(2) (R-R')- p (R') ( 4 ) 

that appears in the vortex average [see Eq. (2)] is a potential determined by the charge 
distribution P(R) for an arbitrary string(s) of elementary dipoles that connect oppositely 
charged probes, p(r) = ±1, of the correlation function (2) (see refs. 13 and 14) . Last, the 
line integral over the vector potential on the right-hand side of Eq. (2) follows the path 
(P) traced out by this track of dipoles. 

Before we go on to compute general phase autocorrelations of the 2D vortex lattice 
using expression (2), it is instructive to first compute the phase rigidity. This quantity is 
known to be given by the one over the dielectric constant of the non-neutral Coulomb gas 
ensemble (3) : 4 ' 15 

p s /J = 1 - lim^o(27r) 2 (Q n Q_ n )/ATk 2 a 2 . (5) 

Here = e Q(R) is the Fourier transform of the vorticity and M denotes the total 
number of sites on the square lattice. Surface barriers are assumed to prohibit "floating" of 
the 2D vortex lattice. 3 ' 4 ' 7 Now suppose that each vortex is displaced by u(R) with respect 
to the zero-temperature triangular vortex lattice. Conservation of vorticity dictates that 
this displacement field is related to fluctuations in the density of the 2D vortex lattice by 

— * 

Q — f = —V • u. This continuity equation is understood to be coarse-grain averaged on a 
scale larger than the triangular vortex lattice constant, a a- Incompressibility is the result: 
rotational invariance gives 

{QkQ-k) = ^ 2 ([Ei^] 2 )/ 2 as k -»■ °- Substitution into Eq. (5) 

then yields the formula 

P£ = 1 _27rJ tt([£'ru] 2 ) 

J k B T ' iV vxa 2 x 1 } 

for the phase rigidity in terms of fluctuations of the center of mass of the vortex lattice, 
where rigid translations are excluded. In obtaining Eq. (6), we have used the trivial identity 
Ma 2 = N vx al xl where iV vx denotes the total number of vortices, and where a vx is the square 
root of the area per vortex. Notice that vortex/anti- vortex excitations not associated with 
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displacements of the zero-temperature vortex lattice are neglected, which shall be assumed 
throughout. Monte Carlo simulation results indicate that this approximation is valid at 
melting. 4 Eq. (6) for the rigidity can be reduced further by expressing the displacement 
field as a superposition of wave and defect components: 7 ' 8 u = u wv + ua{- Waves contribute 
nothing to the degradation of phase coherence, since ^w W v = if rigid translations are 
excluded. We then obtain the final result 

P± = 1 _ ^£ ^(Ej^df] 2 ) m 
J k B T ' iV vxa 2 x 1 > 

for the phase rigidity of a 2D vortex lattice. Notice that the degradation of phase coherence 
is due exclusively to vortex wandering from defects! 

We now resume with the direct calculation of the phase auto-correlation functions 
(2). The first step to take at this stage is to "remove" the net vorticity of the 2D lattice 
through a singular gauge transformation: 5 (e l 52f p '^ = C sw [p]-C vx [p]-e'52R ®° * P l J P ^ dl ^ 
where C vx [p] = (e 52 r S ® $p ) V x represents the vortex component. Here 5Q = Q — Qo is the 
fluctuation in the vorticity with respect to the triangular vortex lattice at zero temperature, 
Qo(R). It is known that C vx [p] is independent of the particular shape of the path(s) P 
connecting probe points p of the correlation function because of the neutrality condition 
J28Q = (see refs. 13 and 14). And since both the correlation function (2) and its spin- 
wave component, C sw , are manifestly path invariant, the phase factor above that remains 
after the singular gauge transformation must therefore also be path invariant. We then 
obtain the formal expression 

( e *I>^ = C sw [p] ■ C vx [p] • e *X>*», (8) 

where the phase <f>o(r) should resemble the zero-temperature configuration. The next 
step is to employ the cummulant expansion, C vx [p] = exp[— ^((52r $Q ' ^p) 2 )vx], for the 
vortex component, which assumes small distortions of the 2D vortex lattice. The latter 
suggests the Taylor expansion J2r 5Q ■ &p = 52'r^' V<E>p for the argument, 9 which yields 
((52 r $Q'&p) 2 )vx = 52a b ^®p\a-(u{a)u(b))-V$>p\b- The cummulant expansion then yields 
the algebraic decay C vx [p] = exp[r7 vx ^^ 2 ^ p(l)ln(ri 2 /ro) p(2)] for the vortex component 
in the asymptotic limit r±2 — > oo and R a b < oo, with a vortex correlation exponent 

' 2 

^ x = 7r ([ "] )/ Ar v X a 2 x , (9) 

R 
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and with the natural ultraviolet scale r ~ a vx . As in the case of the phase rigidity (7), 
the separation of the displacement field into wave and defect components yields the final 
formula 



for the correlation exponent of the vortex component in terms of vortex wandering due to 
lattice defects. Last, Eq. (8) implies that phase correlations in the 2D vortex lattice state 
decay algebraicly, with a net correlation exponent equal to r\ = i] sw + i] vx . Comparison 
with the spin- wave result ?y sw = k B T/2nJ suggests the identity 1] = kBT/2np s between this 
exponent and the phase rigidity. 16 A moment's inspection of Eq. (7) determines that the 
identity indeed holds true in the limit t] v ^ <C ?y sw , when defects are dilute, which is implicit 
throughout. In conclusion, both the present calculation of the correlation exponent and 
the previous calculation of the phase rigidity are in complete agreement. 

To compute the fluctuations in the center of mass due to defects [Eq. (10)], we 
shall now take the continuum limit, / — > 0, and model the 2D vortex lattice (3) as an 
incompressible elastic medium: 10 ' 11 V ■ u = and i? vx = \v J d 2 R{V x u) 2 . Here v rep- 
resents the local 2D shear modulus. Neither vacancy /interstitials nor disclinations are 
important at low temperature. 7 ' 8 We shall therefore compute first the average (|wdf| 2 ) 
due to the presence of a single dislocation pair anywhere inside of the 2D vortex lattice. 6 
Indeed, suppose that a dislocation pair of extent R±2 lies a distance R from the ori- 
gin. It is then easy to show that the displacement there is given asymptotically by 
•Udf(O) = (7ri?) _1 (6 • R)[z ■ (R x -R12)] -R, where ±b are the equal and opposite Burger's 
vectors of each dislocation. We now quote the elastic interaction energy due to such 
a dislocation pair: 6 E 12 = 7r _1 z/6 2 [ln(i? 12 /adf) — |cos20], where 9 denotes the angle 

— * —* 

between R12 and 6, and where aaf represents the core diameter of a dislocation. Af- 
ter averaging |-Udf(0)| 2 over the Boltzmann weight, one obtains the result (|«df| 2 ) — 
(& 2 /4-7r)(?idf (R12) — !^df(-R 2 2 cos 2#))ln(-Ro/adf)- Here n&{ denotes the density of dislo- 
cations, while -Ro represents an infrared cutoff. 17 The logarithmic divergence associated 
with this scale justifies the neglect of higher-order multipole corrections to Wdf(0). It also 
justifies the neglect of the contribution of the autocorrelator (u^ ■ u^ f ) at different points 
to fluctuations (10) in the center of mass. This is due to the fact (i) that — {ud{(a) -Ud{(b)) 
must decay faster than [(|-Udf| 2 )/27rln(i? /ttdf)](«vx/-Ra6) 2 because 1] VX > 0, and to the 




(10) 



R 
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fact (ii) that the former autocorrelator is short range as a result of disordering by the free 
dislocation pairs. Eqs. (7) and (10) then yield a degradation of long-range phase coherence 

jr 1 + - r^^y(^)- (id 

The second term on the right-hand side above is proportional to the renormalization of 
one over the shear modulus obtained first by Kosterlitz and Thouless (KT). 6 The propor- 
tionality constant is set by the estimate vb 2 = (iv/2^/3)J for the shear modulus of the 2D 
vortex lattice, 10 ' 11 with the Burger's vectors limited to basis vectors of the triangular vor- 
tex lattice (b = a a)- We conclude that phase correlations in the 2D vortex lattice reflect 
the same continuous 2D melting transition that is monitored by the macroscopic shear 
modulus! In the limit of a dilute concentration of dislocations, — > 0, this unique 2D 
melting transition occurs at a temperature ksT^ = ub 2 /An, which is about J/14 using 
the previous estimate for v. 10,11 Eq. (11) therefore implies that the correlation exponent 
that characterizes the algebraic decay of long-range phase coherence (2) acquires a small 
value rj m = ksTm^ /2nJ = (287r) _1 just below the melting temperature. 

Standard 2D melting theory also predicts short-range translational order at temper- 
atures just above T m . 6-8 The previous calculation (11) indicates that phase correlations 
are also short range at T > T m . The phase correlation length in the melted state can 
be obtained from the following construction. At low temperature, T < T m , the previous 
algebraic auto-correlation functions (8) can be recovered by taking the appropriate power 
of a partition function describing the system of dislocations, b(R): 

(e'I>*) = <7o + • (Z d m/Zv[0]) J ' /J ■ (12) 

with 18 

r r ' ' i 

Z di [fje\= \V$ exp[-(/e B T/2J / ) ^ (V^ - fje) 2 + 2m ^ • &A»aJ • (13) 

{b(R)} R R 

Above, ift(R) is a real two-component field that ranges over the zero-temperature vortex 
lattice. The vector field if represents tracks of elementary dipoles over the triangular vortex 
lattice that link probe points of the correlation function on the dual honeycomb lattice (see 
refs. 13 and 14), while e represents any unit vector. The exponent J' / J = va\/2it 2 J is 



equal to the ratio of Tm to Tc = ^J/ks, which is approximately equal to (7n)~ 1 by the 
previous estimate for the 2D shear modulus. Last, the pre- factor g ~ 1 is raised to the 
number of positive (or negative) probes, n+, in the correlation function. Residual angular 
interactions between dislocations are suppressed by the vector Coulomb gas (13). 18 It is 
easy to show that the system of dislocations is in a "plasma" phase at high temperatures, 
T > T m . Following the standard prescription, 13 ' 18 the latter indicates that the two-point 
phase auto- correlator (12) decays exponentially above the melting temperature, with a cor- 
relation length £ vx that is proportional to the Debye screening length Ad of the "plasma" . 
The phase correlation length £ vx therefore diverges exponentially as the system of vortices 
begins to freeze exactly like In Ad oc (T — T m )~ v , with exponent v = 2/5 in the case (13) 
when residual angular interactions are suppressed. 18 

The two-point phase correlator is easily calculated directly from the Villain model (1) 
at high temperatures, however. 5 One obtains 

( e i*(l) e -i*(2)) = goe -ri2/^ e ^f; Mr)-df (M) 

asymptotically for this auto-correlation function in the disordered phase, where the phase 
correlation length is inversely related to the line-tension of the strings in the dual theory (1). 
Although this result exhibits exponential decay just like the previous one (12) at T > T m , 
it also exhibits a non-trivial phase factor, unlike the prior result. The resolution of this 
paradox is due to the existence of an intermediate hexatic phase just above the melting 
point, which was discovered by Halperin and Nelson after the initial KT analysis of 2D 
melting. 6 The phase correlations (14) in the vortex liquid phase at high temperatures, 
T > T h: show no vestiges of the zero-temperature vortex lattice asymptotically. The 
intermediate hexatic phase, T m < T <Th, retains long-range orientational order coincident 
with the triangular vortex-lattice state, 19 and this is reflected by the trivial phase factor 
in the phase auto correlator (8). These subtleties are summarized by Table I. 

We shall close by reviewing arguments for the existence of a continuous melting transi- 
tion in the layered 3D XY model with low uniform frustration, which describes the mixed 
phase of layered superconductors. 5 The author has recently performed a partial duality 
analysis of this system in the limit of weak Josephson coupling that leads to the parti- 
tion function Zqg = X!{n z }(^ ' /^1 >2 ) N ^ nz ^iC[pi} in terms of an integer link field n z (r,l) 
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on points r between adjacent layers / and / + 1. Above, C\p{\ represents the phase auto- 
correlation function (2) of an isolated layer / probed at piif) = n z {f, I — 1) — n z (r, /), while 
N[n z ] counts the total number of fluxon charges, n z = ±1. The latter system is dilute 
in the weak-coupling limit reached at large model anisotropy parameters, 7' — > 00. It 
has been shown in this Letter that C[p] is formally identical to the corresponding auto- 
correlator for the 2D XY model without frustration at criticality in the vicinity of the 
2D continuous melting transition (see Table I). Since the layered 3D XY model without 
frustration shows a continuous order /disorder transition, 20 and since this model is also 
described by a partition function like Zcg at weak coupling, we conclude that the layered 
3D XY model with low uniform frustration should show an analogous continuous melting 
transition at weak coupling. 5 

The author thanks Charles Creffield, Angeles Vozmediano, Ed Rezayi and Victor 
Gurarie for discussions. 
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TABLE I. Listed are the phase auto-correlation functions at long-range and the phase 
rigidities obtained here for the various phases of the 2D XY model with a low uniform 
frustration. Rigid translations of the vortex lattice are prohibited. 



Phase 


/ e #(l) e -i0(2)\ 


Ps/J 


Observations 


Vortex Lattice (T < T m ) 
Hexatic (T m < T < Th) 

Vortex Liquid (T > T h ) 


9o (ro/r 12 ) v e^oWe"^ 2 ) 

g Q g-T-12/Cvx e i<^o( 1 ) e -*<M 2 ) 


Vsw/V 





solid 
not rigid 

not rigid 
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